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SUMMARY 


During Che pest year, Che granC aupporCed research on several 
aspecCs of crusCal deforaaclon. The relaCion beCween earChquake dis- 
placemenCs and faulc dimensions was sCudied in an efforc Co find 
scaling laws chaC relaCe sCaCic parameCers such as slip and sCress 
drop Co Che dimensions of Che rupCure. Several implicacions of Che 
sCaCic relacions for Che dynamic properCies of earChquakes such as 
rupCure velocicy and dynamic sCress drop were proposed. 

A Cheorecical basis for earChquake-relaCed phenomena associaCed 
wich slow rupCure growch or propagaCion, such as delayed mulciple 
evencs, was developed using Che sCress inCensiCy facCor defined in 
fracCure mechanics and experimenCal evidence from sCudies of crack 
growCh by sCress corrosion. 

Finally, exCensive sCudies by Japanese geologisCs have esCab- 
lished Che offseC across numerous faulCs in Japan over Che lasC one 
hundred Chousand years. These observaCions of inCraplaCe faulcing are 
being used Co escablish Che spaCial variacions of Che average v 'rain 
race of subregions in souChern Japan. HisCoric inCraplaCe seismiciCy 
has also been inCerpreCed as a regional sCrain paccern for comparison 
wich Che long-Cerm geologic paCCern. The hisCoric seismic daCa were 
used Co esCablish a relaCioii becween Che lengCh of a faulc and Che 
maximum size earChquake possible on Che faulc. The geologically 
observed slip race Chen was combined wich Che observed size disCribu- 
Cion of earChquakes (b-value) Co predicC repeaC Cimes beCween inCer- 
vals of differenC inCensiCy ground shaking, or one form of seismic 
risk map. 
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COMPLETED RESEARCH 


Scaling Laws for Large Earthquakes 

It is observed that i.'he mean slip in large earthquakes correlates 
linear with fault length L and is not simply related to fault width, 
W. If we interpret this in terms of an elastic model, it implies that 
static stress drop increases with aspect ratio (L/W). He also observe 
a tendency, both for strike-slip and thrust earthquakes, for aspect 
ratio, and hence static stress drop, to increase with seismic moment. 
Dynamic models of rupture of a rectangular fault in an elastic medium 
show that the final slip should be controlled by the fault width and 
scale with the dynamic stress drop. The only way these models can be 
reconciled with the observations is if dynamic stress drop correlates 
with fault length so that it is also nearly proportional to aspect 
ratio. This could only happen if fault length is determined by the 
dynamic stress drop. There are several serious objections to this, 
which lead us to suspect that these models may be poor representations 
of large earthquakes. Firstly, it conflicts with the observations for 
small earthquakes (modeled as circular sources) that stress drop is 
nearly constant and independent of source radius. Secondly, it con- 
flicts with the observation that fault length is often determined by 
the size of seismic gap in which the earthquake occurred. We specu- 
late that the boundary condition at the base of the fault, that slip 
is zero, is unrealistic because that edge is in a ductile region at 
the base of the seismogenic layer. In a madel in which slip is not so 
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constrained at the base of the fault nor at the top (the free sur- 
face), such that no healing wave originates from these edges, final 
slip would be determined by fault length. The observations would then 
be interpreted as meaning that the static and dynamic stress drops of 
large earthquakes are nearly constant. Geodecic data for the 1906 San 
Francisco earthquake show that this model iu not satisfactory, 
either. These two extreme models predict very different scaling of 
the dynamics of large earthquakes. . The width-dependent model predicts 
that average particle velocities are larger for long ruptures but the 
rise time will be the same as in a shorter event of the same width. 
The length-dependent model predicts the opposite. 

Time-Dependent Rupture 

Using the concepts of fracture mechanics, we develop a theory of 
the earthquake mechanism which includes the phenomenon of suberitical 
crack growth. The theory specifically predicts the following phe- 
nomena: slow earthquakes, multiple events, delayed multiple events 

(doublets), postseismic rupture growth and afterslip, foreshocks, and 
aftershocks. The theory also predicts that there must be a nucleation 
stage prior to an earthquake, and suggests a physical mechanism by 
which one earthquake may 'trigger' another. 

These predictions are obtained by combining two fundamental con- 
cepts. The first is that 


Ic ■ cAtiGc 
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and the second, that 


k 


„ ,X vl/n 
o 


where k is thw stress intensity factor, 4 t stress drop, X rupture 
length, X rupture velocity, C a geometrical factor, and Kq, Vq, 
and n are material constants. The first is a fundamental result of 
fracture mechanics; the second describes stress corrosion cracking, a 
well-established physical process. 

We investigate in detail two phenomena of special interest and 
which are not predicted by ordinary fracture mechanics: nucleation 

and delayed multiple events. In the first case, we find that all 
earthquakes must be preceded by quasistatic slip over a portion of 
their rupture surfaces, but it may be difficult to detect in prac- 
tice. In the second case, we studied two pairs of delayed nultiple 
events that were separated by the same 'barrier' in order to calculate 
Kq and n. The results are that the stress corrosion index, n ■* 24 
and the stress corrosion limit, 10^ < Kq < 10^ bars km^/2. 
The latter uncertainty is due to uncertainty in Vq, which we assume 
to be in the range 10“^^ Vg < 10”^ km/hr. 
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Japanese Intraplate Crustal Deformation 

The his tor ic^*! record of large (M _> 6.9) earthquakes anrt geo- 
logically determined rates of slip on Quaternary faults in intraplate 
Japan (Honshu and Shikoku) are used to estimate the average rate of 
seismic moment-release (Mj^j) for the last 400 years and during the 

a, 

late Quaternary, respectively. Values of estimated from the two 
data, sets are nearly equivalent in regions where seismic activity is 

e_ 

concentrated on land. This observation suggests that in intra- 

plate Japan has been constant during the late Quaternary and is rel- 
atively free from secular variation when averaged over periods of sev- 
eral hundreds of years. The fabric of Quaternary faults in northeast 
Japan indicates that the principal axis of horizontal compression is 
oriented eastward, parallel to the relative plate velocity between the 
Eurasian and Pacific plates along the Japan trench. Similarly, the 
pattern of Quaternary faulting in Shikoku shows a compressive stress 
field that trends northwest approximately parallel to the relative 
plate velocity vet'ttor between the Eurasian and Philippine Sea plates 
along the Mankai trough. The mode of intraplate deformation in 

Shikoku is, however, in marked contrast to that observed in northeast 

s 

Honshu. Virtually all in Shikoku may be attributed to slip 

s 

along one fault, the Median Tectonic Line (MTL) . In contrast, M£j 
in northeast Honshu is divided among a regional set of reverse-type 
faults that strike northerly, perpendicular to the direction of plate 
subduction along the Japan trench. The easterly strike of the MIL is 
consistent with the idea that oblique plate subduction along the 
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Nankai trough is partially accommodated by right-lateral motion along 
the MIL. Conversion of Mij in northeast Japan to strain-rates 

a 

(c£j) suggests that about 3 to 9X of the relative plate motion 
between Japan and the Pacific plate (" 9,7 cm/year) is accomoodated as 
a permanent eastward shortening (3-9 rao/yr) of intraplate Japan. 
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Figure I . Regional averages of the rates of moment-release, hori- 
zontal compressive strain and crustal shortening calculated from the 
40D year record of seismicity, Principal directions of horizontal 
compressive strain and shortening are denoted by dashed bars. Note 
that the values in the Izu region are calculated from seismicity 
within the small dashed box that encloses the Izu Peninsula. 

Figure 2. Regional and grid-by-grid atverages of the rates of moment- 
release, horizontal compressive strain, and crustal shortening calcu- 
lated from geologically determined slip-rates of Quaternary faults. 
Principal directions of horizontal compressive strain for the larger 
intraplate regions and smaller grids arc shown by dashed bars and thin, 
solid lines, respectively. Smaller numerals within the region of cen- 
tral Japan give average rates for the local areas enclosed by the 


dashed boxes. 
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CONTINUING RESEARCH 
Gruatal Deformation in Alaska 


The grant has partially supported transmission via satellite o£ 
geodetic data from two seismic gaps in Alaska^ It has also aided the 
development of pressure recording systems used to detect long-term 
water level changes resulting from tectonic displacement of the 
earth. Seven 17 m long carbon-fiber strainmeters were installed in 
shallow surface sites near the coast of the Gulf of Alaska within the 
Yakataga seismic gap. All the instruments are located above the fore- 
cast rupture zone. Hif<arly-mean strain data are telemetered via the 
Landsat 3 satellite for analysis in New York. 

Transmissions have been received from six of the strainmeters 
since November 1980. One transmitter and two strainmeters h&d early 
malfunctions and a third strainmeter ceased operation in January 
1981. Three strainmeters continue to transmit data of good quality. 
The least-count resolution is approximately 10~^ strain. 

The load tide signal has a maximum peak amplitude of 4 x 10~^ 
strain. Superimposed on the tidal signal are strain rates that vary 
from 10'~^/day lasting for as long as 30 days and strain rates of 5 
X 10~^/day lasting less than ti days. The higher strain rates are 
probably associated with temperature, rainfall, and ground expansion 
during freezing processes. Our ability to distinguish tectonic sig- 
nals amid these large atmospherically-induced strains depends on the 
magnitude and spatial geometry of likely precursory strainfields. 
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since itoim conditions ere essentially synchronous across the Yakataga 
coastline, ve cannot use temporal coherence alone to search for possi- 
ble tectonic sis‘;als within the strain time series. In order that we 
may detect unambiguously a strain precursor, the magnitude of the pre- 
cursory tectonic strain rate nust either exceed the maximum observed 
strain-noise rate (e.g., l0”^/day) or we naisf be able to dio*"if fwish 
the precursory strain signal from the noise by virtue of it > ^ 

geometry. The nature of a short-term strain precursor to a magnitude 
8. 1-8 ,3 earthquake cannot be surmised from historical data since no 
continuous strain data from within the epicentral region of an event 
of this magnitude have been reported. Moreover, the spatial form of 
strainfields generated by fast-moving storms are likely to be diffi- 
cult to predict. Because of the large size of strain precursors that 
have been reported for twenty out of forty documented precursory 
strain events (Bilham, 1981), we remain optimistic about our detection 
capability. 

A newly designed sea-level instrument is to be installed in 
Alaskan seismic gaps in 1981 to monitor crustal movements associated 
with future seismicity. The device consists of two ceramic pressure 
gauges: one to monitor barometric variations and the other installed 
in shallow water on the coast to measure sea-level variations. The 
instruments have a resolution of approximately I mm over a range of 
more than 10 m. Sea-level, barometric pressure, and sea water temper- 
atures will be averaged for 12 minutes, stored for up to three hours, 
and telemetered via the GOES satellite. In each region we will 
install two permanent gauges and one portable gauge. The portable 
gauge will be operated for variable lengths of time in order to estab- 
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lish several vertical-control points within and near each seismic 
gap. These control points will be reoccupied yearly. 


Cruatal Deformation in Iceland 

The South Iceland Seismic Zone is an east-west trending zone of 
large historic earthquakes that transfers motion from the submarine 
Reykjanea spreading center to the Eastern Volcanic Zone in southern 
Iceland. Major earthquake sequences affecting most of the 70 km Inng 
zone recur at intervals of 50 to 100 years. Earthquakes of magnitude 
7 or greater have been recorded. No major earthquake sequence has 
occurred since 1S96, and the whole zone has had a low level of seismic 
activity instrumentally located in it for the past 50 years. Thus 
there is a high probability for the occurrence of a large earthquake 
in the zone in the next few decades. The grant has helped to initiate 
a program for installing 5 water level sensors in Lake Hestvatn, 
located in the center of the seismic zone. The lake surface is used 
as a reference level to measure crustal tilt in the area. Anomalously 
rapid tilting has been observed before earthquakes in some cases. 
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